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Abstract. We discuss the construction of cosmological models within the framework 
of Macroscopic Gravity (MG), which is a theory that models the effects of averaging 
the geometry of space-time on large scales. We find new exact spatially homogeneous 
and isotropic FLRW solutions to the MG field equations, and investigate large-scale 
perturbations around them. We find that any inhomogeneous perturbations to the 
averaged geometry are severely restricted, but that possible anisotropics in the corre- 
lation tensor can have dramatic consequences for the measurement of distances. These 
calculations are a first step within the MG approach toward developing averaged cos- 
mological models to a point where they can be used to interpret real cosmological data, 
and hence to provide a working alternative to the "concordance" ACDM model. 
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Introduction 

The currently favoured "concordance model" of cosmology treats the entire observ- 
able universe as a single Friedmann-Lemaitre-Robertson- Walker (FLRW) solution of 
Einstein's equations. The complex hierarchy of structures that astronomers observe 
is then accommodated by allowing for small fluctuations to this geometry, and a lin- 
earized version of Einstein's equations is employed to determine their behaviour. The 
mathematical simplicity of this approach allows extremely complicated distributions of 
matter to be considered, and permits one to perform detailed calculations of a variety 
of complicated physical phenomena. It relies, however, on the assumption that the 
equations that govern the large-scale 'average' geometry of the Universe are the same 
as Einstein's equations. This will not be true in general. 

What we have in reality is a complicated fitting problem: We make observations 
over various different scales (depending on the particular observable in question), and 
we then look for the FLRW model that best fits either one or all of these observables. 
Ideally it would be better to make direct observations of the matter distribution, and 
then solve Einstein's equations with some reasonable set of boundary conditions. Un- 
fortunately this is not currently possible, due to the complexity of the equations, and 
so we are forced to revert to the former approach. In order for this procedure to have 
any meaning, however, we must be prepared to answer a series of questions: What 
is the relation (if any) between the fitted FLRW model and the actual geometry of 
space-time? What are the evolution equations that the fitted FLRW model should 
obey? Should the FLRW models we fit to observations made on different scales be 
expected to be the same, and, if not, how (if at all) are they related? 

One way to go about trying to answer these questions is to attempt to develop and 
understand procedures for averaging the geometry of space-time, thereby providing an 
explicit link between the microscopic geometry of space-time, and the macroscopic 
averages that are used in cosmology. This is the goal of Zalaletdinov's theory of 
Macroscopic Gravity (MG) [1], which consists of a prescription for averaging tensors, 
vector and scalars in a covariant way, as well as providing a set of field equations that 
the averaged geometric quantities are expected to obey. Despite the difficulties in MG 
[2-5, 12], it does serve as a framework within which calculations can be performed. 
Highly symmetric exact solutions to the field equations of MG have been found in [2, 
3, 6, 7], the prescribed averaging procedure has been directly applied to space-times 
that are already close to being spatially homogeneous and istropic in [8-10], and a 
recent attempt to use the exact solutions found in [3] and [6] to interpret cosmological 
observations has been made in [11]. For other approaches to averaging in cosmology 
the reader is referred to [12, 13], and references therein. The relationship between MG 
and some of these approaches is considered in [8]. 

Our aim here is to further develop the construction of cosmological models within 
the theory of MG, thereby building on previous work that has identified exact spatially 
homogeneous and isotropic solutions to the MG field equations. In particular, we aim 
to investigate cosmological perturbation theory within MG theory, as well as making 
some first steps towards understanding the propagation of rays of light. In Section 
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1 we briefly summarize the essentials of MG theory, and then proceed to discuss its 
FLRW solutions in Section 2. We first recap some known solutions, before presenting 
new solutions that have not previously appeared in the literature (details of these new 
solutions can be found in the appendix). In Section 3 we then investigate perturbations 
around the known FLRW solutions, finding that MG theory requires that either: (i) 
The background macroscopic geometry must be the same as the FLRW solutions of 
Einstein's equations (up to the order of perturbations considered), or (ii) The pertur- 
bations themselves must be spatially homogeneous (or take a very restricted form). We 
then proceed in Section 4 to calculate measures of distance in the averaged geometry 
by using the average of the Ricci and Weyl curvature tensors in the Sachs equations. 
We find that in the simplest solutions the correlation tensor acts as a spatial curvature 
term in the Sachs equations, in the same way that it does in the MG field equations. 
However, this result does not extend to more general FLRW solutions which exhibit 
considerable freedom. Finally, we conclude with some brief remarks. 

1 Macroscopic Gravity Theory 

Macroscopic Gravity (MG) is Zalaletdinov's non-perturbative approach to describing 
the behaviour of space-time on large scales [1]. We will briefly summarize it here. 

1.1 Covariant Averaging 

Zalaletdinov's approach involves averaging the geometric objects that exist on the 
space-time manifold, and constructing field equations for these averaged quantities 
based on averaging the Cartan equations for the microscopic geometry. The first step 
in this is to provide an averaging procedure that is covariant, and that maintains the 
tensorial properties of the object that is being averaged. To do this, Zalaletdinov 
defines that the average of an object p a p , over some closed region of space-time £ 
that contains the supporting point x, to be [1] 

(/'•;:::(•'•)> = ±- I ^d*xyJr(x>)A%,(x, x')a»' p (x, x>)..., (1.1) 

where Vs = J s ^J—gd^x is the 4- volume of S, and A a fM r(x, x') and A v Jx, x') are bilocal 
averaging operators. These operators can be written as the product of vector bases at 
two different points in E as A a p(x,x') = e a i (x')e l /3 (x), where the structure functions, 
Cij k of [e-i,ej] = Cij k ek (i.e. the coefficients of anholonomity) , are constants. 

With this definition we can now consider the average of various geometric objects. 
A key point to bear in mind here is that averaging non-linear quantities is not, in 
general, the same thing as constructing those same quantities from the average of their 
arguments. One must therefore proceed with care. Let us specify that we will write the 
un-averaged metric, connection and Riemann tensor in the normal way as g pu , r M and 
R^ upa - We can then write the averages of these quantities as (g^), (T^ ) and {R^ upa )- 
It is now, for example, no longer the case that the curvature tensor constructed from 
(T^ up ) is given by ). Instead, following Zalaletdinov, we denote this object as 

M» ua e = 0«<T%) - dp{T» va ) + (T\ a ){T\») - (T»)(T% a ), (1.2) 
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where M ua/3 ^ (R va p), in general. Likewise, it is assumed that there exists a non- 
metric compatible and symmetric connection, IP [1], such that 

(R^) = d a ii% - d^w va + ip^if^ - n%rr ra , (i.3) 

where IP 7^ (r^p)) * n general. The difference between these curvature tensors, 

(Tuafi = ^uj ~ (1.4) 

^(rV^l^rj^,), (1.5) 

where under-lined indices are not included in anti-symmetrization, is known as the "po- 
larization tensor" . The difference between the averaged connection and the connection 
of the averaged Riemann tensor yields an "affine deformation tensor" : 

^Hrg-nv (1.6) 

In what follows we will often use an over-bar on a quantity to denote the average, 
rather than angular brackets, as this will allow us to be more concise. 

1.2 The Macroscopic Field Equations 

Under the assumption that the connection = (r M ) is compatible with the metric 
9iw = {9^u)i an d that we assume the following splitting rules for products of connection 
and metric, (T» up g%) = (T» vp ){g%) and (T^T^g^) = (T a ^T^) (g%) , one can 
then construct the field equations of MG: 

/ e M 7/3 - U%r v M, v = SttGT^ - (Z%, 7 - U%Q, v )r\ (1.7) 

where T € 7 is the averaged energy-momentum tensor, and where Z a ^ = 2Z a e v p and 
Q^ u = Z a ^ ua . Here Z a ^ ua is known as the "connection correlation tensor", and is 
given by 

We note that the "polarization tensor" is the trace of the "connection correlation 
tensor", via Q a p pu = — 2Z e « We assume that = (T^ u ) can be written as a 
perfect fluid, such that 

n=(p + p)u\+#, + ^ I (1-9) 

where p, p and Y^ v are the energy density, isotropic pressure, and anisotropic stress of 
the fluid, respectively. The 4-vector ?p exists in the manifold with averaged geometry 
g~n V , and is the velocity 4-vector of the fluid. 
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1.3 The Connection Correlation Tensor 

From the definition of the connection correlation tensor (1.8), it can be seen that 
Z a p£ va ) = 0, Z a ^ l va = —Z^"^, and that we have the cyclic constraint 

The connection correlation tensor is also assumed to satisfy the "equi-affinity" con- 



straint 

and the differential constraint 



Z\^ va = (1.11) 

where 1 1 denotes a covariant derivative with respect to the averaged connection, T tJ- up . 
The integrability conditions for this equation are 

Z € /3[fi' y 6vM a £ K7T ] — Z a e [^SuM e p K7r ] + Z a p[ l j, e syM y e / v r ] — Z a ^ ' & M € $ K7T ] = 0. (1-13) 

With the assumption of vanishing 3 and 4-form correlation functions for the connection 
and assuming Eq. (1.12), we have the quadratic constraint 

ZS 6 rycx I ry8 fj, ryQ a 

/3[7 ktt Se va\ ~r £ /3[ 7 vp£ kit Se] 

I ryOL S ryp , ryOL fj, ryO S 

~r 6 /3[7 vrj^ Se ktt] T & /3[j 5e^ ktt up] 

+ Z a pyf SeZ^ vp 5 ktt\ + Z a p] a S j ^Z s/'va] — 0- (1-14) 
Eqs. (1.12)-(1.14) are then solved to determine the connection correlation tensor. 

1.4 The AfRne Deformation Tensor 

The affine deformation tensor must obey the differential constraint 

1 
2 

From the Bianchi identities, and Eqs. (1.4) and (1.12), we then have R a m pr7 \\x\ — 0, 
where R a p pa = {R a i3 pa ), which gives 

A^ [p R\ aX] -A a e[p R^ aX] =0. (1.16) 

Eqs. (1.15)-(1.16) are solved to determine the affine deformation tensor. 
2 Macroscopic FLRW Solutions 

The FLRW solutions to Eqs. (1.7)-(1.16) have been studied by Coley, Pelavas and 
Zalaletdinov [6], and van den Hoogen [3]. The latter of these gave a systematic and 
explicit investigation of all spatially flat solutions in which g pu , Z a ^ urj) and A^ va are 
invariant under the full six-dimensional group of Killing vectors that characterize the 
spatially homogeneous and isotropic FLRW geometries. We will now briefly outline 
these solutions (which we will call Type I), as well as some more general FLRW so- 
lutions in which the connection correlation and affine deformation tensors no longer 
have the same symmetries as g pu (which we will call Type II). These latter solutions 
have not appeared in the literature before now. 
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2.1 Type I Solutions 

These solutions have correlation and affine deformation tensors that are invariant under 
the same group of Killing vectors as g^ u . In their most general form they are fully 
specified by three independent constants: A, h 2 and b\. All three of these parameters 
appear in the (constant) components of the connection correlation tensor, with the 
affine deformation tensor being fully specified by A only (see [3] for details). 
Writing the macroscopic line-element as 

ds 2 = g^dx^dx" = -dt 2 + a 2 (t) [dx 2 + dy 2 + dz 2 ] , (2.1) 

we then have that the macroscopic field equations (1.7) give 

er 3 a z 

where over-dots denote differentiation with respect to t, and where p is the macroscopic 
energy density from Eq. (1.9) that obeys 

p + 3^(p + p) = 0. (2.3) 

The macroscopic geometry therefore evolves like the FLRW solutions of Einstein's 
equations with an extra spatial curvature term in the macroscopic equivalent of the 
Friedmann equation, even though the macroscopic geometry is spatially flat. 

The constants h 2 and by do not affect the dynamics of the averaged space-time, 
and with h 2 = bi = the non-trivial components of the connection correlation and 
affine deformation tensors are 

2 

yj k yj ft yj j yk k _ _ A 2 Icy A\ 

33 jk kj kk ^ jj kk ^ jj kk 10 / 

yi k yj i yi i yj k _ j\2 /n r\ 

33 ik ^ ij kk ^ jj kk ^ ij ik 10 \ / 

where [ijk] is chosen to be one of the ordered triples {[1, 2, 3], [2, 3, 1], [3, 1, 2]}, and 

A\ l = Aa{t), A i i0 = A° 00 = ^- J (2.6) 

a[t) 

where i takes on values 1, 2, 3. The observational consequences of a simple application 
of this solution have been investigated by Clarkson et al. [11]. 

2.2 Type II Solutions 

Further solutions with a FLRW macroscopic metric, g^, can be found if Z a '^ va and 
A^ ucr are allowed to be less symmetric. From Eq. (1.7) it can be seen that this is 
possible because only certain contractions of these tensors are required to admit the 
full six- dimensional group of Killing vectors that are required for g^ v to be an FLRW 
metric. That is, the correlation and affine deformation tensors do not have to be exhibit 
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invariance under all spatial translations and rotations themselves in order for g^ v to 
be FLRW. 

We find that, by relaxing the symmetry requirements on Z a ^ va and A^ V(T , there 
exists a family of four solutions to the algebraic constraint Eqs. (1.13), (1.14) and (1.16) 
that are parameterized by the ordered pair (a 1; a 2 ) G {(1/2, 0), (1/2, 1/2), (1, 0), (1, 1/2)}. 
In this case it is found that Z a ^ UCT contains: 

51 functions fi(x,y,z), 2 functions J-i(x,y), and 1 constant A, (2.7) 

and A^ va . contains: 

2 functions £>j(x, y) and 2 constants A, B. (2.8) 

where A is the same constant as found in the Type I solutions, in the appropriate 
limit. Equation (1.12) then provides a lengthy set of differential equations that must 
be obeyed by the fi(x,y,z), and Eq. (1.15) provides the differential constraints (5.1), 
which essentially define the connection correlation functions Fi{x,y) in terms of the 
functions Bi(x,y) found in the Affine Deformation tensor. 

The macroscopic field equations (1.7) then reduce to the same form as in the Type 
I solutions described above, such that Eqs. (2.2) and (2.3) must be satisfied. That 
is, the only quantity that enters into the macroscopic metric g^ v from the correlation 
tensor is the single constant A, which appears as a spatial curvature term (this can 
be seen to be due to assumption 4 in the appendix). The extra freedom introduced 
into the macroscopic description of space-time by allowing Z a ^ u(j and A^ ua to be less 
symmetric than does not therefore, in this case, result in any different behaviour in 
the macroscopic metric. It does, however, result in extra terms in both the connection 
correlation tensor and affine deformation tensor that we will find can be important for 
observations made within the space-time. 

The result that the macroscopic metric obeys Eq. (2.2) is true independent of 
whether or not we set the fi(x,y,z) to be constants. For further details of this gener- 
alized solution, the reader is referred to the appendix. 

3 Perturbed Spatially Homogeneous and Isotropic Solutions 

We now want to look for solutions to the equations of macroscopic gravity, as outlined 
in Section 1, that are less symmetric than the spatially homogeneous and isotropic 
macroscopic geometries that were studied in Section 2. We will do this by allowing 
for small fluctuations about the FLRW geometries that we already know, and then 
performing a perturbative expansion in the fluctuations at the level of the equations 
that govern g^ u , Z a ^ V(T and A^ va . This method is routinely used when studying the 
'almost' FLRW solutions that are heavily relied upon in the standard cosmo logical 
model. It is within this type of frame-work that predictions involving weak lensing, 
galaxy correlation functions, and CMB observations are routinely calculated. As the 
MG equations describe macroscopic behaviour, we shall only consider fluctuations on 
scales comparable to or larger than the averaging domains E. In what follows we will 
refer to these as 'macroscopic scales'. 
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3.1 Perturbed Macroscopic Variables 

Let us begin by writing the perturbed macroscopic geometry as 

ds 2 = a 2 (r]) [- (1 + 2$) drj 1 + (1 - 2*) (dx 2 + dy 2 + dz 2 )] , (3.1) 

where rj = J dt/a{t) is the conformal time, and is used here to simplify the equations. 
The perturbed velocity 4-vector can then be written as 

= -(l-<S>,v i ) , (3.2) 

a v ' 

where the appearance of $ in the temporal component is due to the normalization 

u^u^ = —1, and v l is the peculiar velocity. For the perturbed energy- momentum 
tensor we then write 

T° = -Sp (3.3) 

T*i = -{p + p)Vi (3.4) 

T) = S t J Sp+(p + p)D i j E, (3.5) 

where = didj — SijV 2 is a traceless spatial derivative operator, Sp is the energy 
density perturbation, 5p is the pressure perturbation, and £ is the scalar anisotropic 
stress. We further assume that we can write 

= ^ (0) V- + 5Za M - ( 3 - 6 ) 
= A^ va + 5A\ a , (3.7) 

where a superscript (0) denotes the value of a quantity in the exact FLRW solution that 
is being perturbed, and 5Z a ^ ^ va and SA 1 *^ indicate small fluctuations. In general, 

^Z a ^ va has 720 independent components, and 5A^ ua . has 40 independent components. 

3.2 Perturbed Equations 

We now want to substitute our perturbed quantities (as specified in Eqs. (3.1)-(3.7)) 
into Eqs. (1.12)-(1.16), and to expand the result in some 'order of smallness' parameter, 
e. In doing this, we take our fluctuations to be of the following order: 

$^^^^^5p^^^E~ 5Z a p/ ua ~ SA"^ ~ e. (3.8) 

The O(0) part of these equations is automatically satisfied, if we perturb around the 
exact FLRW solutions that we presented in Section 2. We are then left with a set of 
equations that the fluctuations must obey. 

The general form of the perturbed equations is very lengthy, with a large number 
of functions describing the fluctuations (in particular, in bZ a ^ va and bA^ va ). The 
variables that we are most interested in, however, are those that are involved in the 
macroscopic geometry and the macroscopic fluid description of the matter, i.e. $, 
Sp, Sp, £ and v % . Information can be gained on the form that these fluctuations must 
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take in the general case, but in order to illustrate the form of the perturbed equations 
in this paper we will first consider some simplifying special cases. 

Let us begin by considering small fluctuations around the Type I solutions from 
Section 2.1. We can set some of the fluctuation in the correlation tensor to zero by 
choosing 

In the language of [3], this corresponds to setting the electric components of the corre- 
lation tensor to zero. This condition means that the quadratic constraint, Eq. (1.14), 
is automatically satisfied. It also reduces the number of independent components in 
5Z° 1 p^ va . to 360. The cyclic and equi-affine constraints, Z 01 ^ , = and Z a a ^ va = 0, 
together with the integrability conditions, Eq. (1.13), and the differential constraint 
equation of the affine deformation tensor, Eq. (1.15), then leave only 64 independent 
variables in 5Z a ^ U(j . 

3.2.1 The Macroscopic Field Equations 

After lengthy calculation, the macroscopic field equations (1.7), at first order in per- 
turbations, are found to give 



7^-^ - 7^? ) [$ - * - 87iGa 2 (p + p)E] = 



(3.10) 



[$ - * - 8vrGa 2 (p + p)E] = 0, (3.11) 



dx 1 dxi 

where % ^ j, and no sum is implied in Eq. (3.10). We also find 

V 2 $ - m ($' + W$>) = AnGa 2 (5p + 5p grav ) (3. 12) 

A ($' + m) = AixGa 2 (p+p- Vi (3.13) 

$" + WS' + 2U& + {2W + H 2 ) * + ^V 2 (^ - $) = AnGa 2 Up - ^S Pgr J\ , (3.14) 

where a prime denotes differentiation with respect to rj, and we have defined % = a' /a. 
In Eqs. (3.12)-(3.14) we have introduced the new quantity Sp grav , which contains all 
the first-order contributions from the perturbed correlation tensor, with the exception 
of the term containing A in Eq. (3.13). 

3.2.2 The Connection Correlation Tensor 

Moving on from the macroscopic field equations, the integrability condition (1.13) also 
means that <3> must obey the following equations: 



A 2 



<9 2 $ <9 2 $ 



[no sum implied] (3.15) 



X W = (316) 
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where % 7^ j and i, j G {1, 2, 3}. The differential constraint equations for the correlation 
tensor, (1.12), then give 



A 2 



dvi dvj 



[no sum implied] (3-17) 



A* = 0, (3,8) 

where i 7^ j. Note that we have not yet assigned an order of smallness to A, so if 
A ~ 0(e 1 ^ 2 ) or smaller then Eqs. (3.15)-(3.18) are automatically satisfied. 

Beyond these equations, we also gain from Eqs. (1.12) and (1.13) constraint 
equations for ty' and v[, as well 42 linearly independent equations involving spatial 
derivatives of the first-order contributions to the correlation tensor. 

3.2.3 The Affine Deformation Tensor 

At first order in perturbations there are initially 40 independent components to 8A a ^. 
The linearized version of Eq. (1.16), assuming Eqs. (3.15) and (3.18), fixes all but 12. 

3.3 Macroscopic Perturbed FLRW Solutions 

Having summarized the equations for first-order perturbations about an FLRW back- 
ground in MG, let us now turn to their solutions. 

3.3.1 The Macroscopic Field Equations 

The off-diagonal spatial components of the macroscopic field equations, given in Eqs. 
(3.10) and (3.11), can be seen to take exactly the same form as they do in the linearly 
perturbed FLRW solutions of Einstein's equations. These equations are usually taken 
to imply that $ — ^ = 87rGa 2 (p + p)£, so that when the anisotropic stress vanishes we 
have $ = ^. In the most general situation, we can see that $ — ^ — 8nGa 2 (p + p)S 
must be of the form 

ci + c 2 x + c 3 y + c A z + c 5 (x 2 + y 2 + z 2 ), (3.19) 

where C\, c 2 , C3, C4 and C5 are all functions of 77 only, and are of first order in our 
perturbative expansion. 

Turning to Eqs. (3.12)-(3.14) we again see that these equations take a similar form 
to the familiar ones that are arrived at when considering perturbed FLRW solutions 
of Einstein's equations. The contributions to these equations from the correlation 
tensor occurs solely in the appearance of the terms containing Sp gvav and A. This is a 
considerable increase in simplicity over the large number of variables that are required 
in order to specify the correlation tensor itself. Moreover, it can be seen that the 
contribution of the correlation tensor to the macroscopic field equations takes the form 
of a perfect fluid with effective energy density and pressure of the form 

3^42 

Pgrav = 5- + fygrav + 0(e 2 ) (3.20) 

A 2 1 

Pgrav = — - o<Wav + 0(e 2 ). (3.21) 
Ct o 



- 10 - 



Clearly, this effective ffuid has the following equation of state: 

Pgrav = -gPgrav + 0(e 2 ). (3.22) 

This is same effective equation of state as that of a spatial curvature term in the 
Friedmann equations, and so the result that the correlation tensor mimicks spatial 
curvature is found here to be valid not only at the level of the background, but also to 
at least first order in perturbations around the background. This extends the findings 
of [6] and [3] beyond perfectly spatially homogeneous and isotropic solutions. 

3.3.2 The Connection Correlation Tensor 

The differential constraint (1.12), and its integrability condition (1.13), provide us at 
first order in perturbations with Eqs. (3. 15)- (3. 18). If we write V\ = + v\ , where 
is the curl- free part of Vi and v!> is the divergence- free part of Vi, then we can 
also write = dG/dx*. If A ~ 0(1) then Eqs. (3.15)-(3.18) show that 6 and <£> 
must take the form given in Eq. (3.19), and if $ = ^> then the same must also be 
true for These equations are therefore very constraining on the form that \P and 
the scalar part of Vj, can take. In fact, if we were to perform the usual treatment of 
these quantities in terms of Fourier modes, then the resulting equations would suggest 
that all modes vanish to first order in perturbations. That is, they would show that 
no inhomogeneous perturbations are possible at all. 

3.4 Interpretation 

First of all, one may consider what happens to the perturbations when the assumptions 
we have imposed on the solution are relaxed. For example, we could have used the more 
general Type II solutions found in Section 2.2. Alternatively, we could have attempted 
to relax our assumptions involving the vanishing of the electric part of Z, as prescribed 
by Eq. (3.9). It is our contention that in these cases one still obtains equations similar 
to (3.15)-(3.18). This is due to the form of Eq. (1.13), which consists of terms built 
from contractions of the correlation tensor, Z, with the Riemann tensor constructed 
from the averaged connection, M. At first-order in perturbations we will therefore 
generically have equations containing terms that contain the zero-order part of Z and 
the first-order part of M. It is exactly these terms that lead to Eqs. (3.15) and (3.16), 
and which should therefore be considered generic. We note here that an exceptional 
case is the one in which the background space-time is Minkwoski space, as in this case 
the zero-order part of Z vanishes, and Eqs. (3. 15)- (3. 18) are satisfied identically. This 
is the reason that it was possible to find non-trivial spherically symmetric solutions to 
the MG equations in [2] (see also [7]). 

It could perhaps be argued that Eqs. (3.15)-(3.18) only need to be satisfied 
in some 'averaged' sense, but this is incompatible with the approach taken here, and 
would, presumably, require some adjustment of the macroscopic theory. The equations 
we have obtained consequently must be taken as constraints on the form of the macro 
metric (3.1), and as such need to be interpreted. 
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If we assume Z can be perturbed as in Eq. (3.6), and that the MG approach is 
generally valid, then we see 2 general possibilities: 

1) . A ~ 0. That is, if Z^ arises from the zero-order solution for the macroscopic 

metric in Eq. (3.1), then the differential constraints for Z given in Eq. (1.12) 
imply that the contribution from the correlation to the zero-order part of the 
macroscopic field equations must be 0(e 1 ^ 2 ) or smaller. Either this, or there 
must exist as yet unknown solutions of the macroscopic field equations for Z^ 
(that is, different to the form of Z^ given in Section 2, which may be possible 
but not proven). 

2) . $ = \P = 0. That is, the macroscopic metric in Eq. (3.1) can only contain spa- 

tially homogeneous perturbations (or, at most, the very highly restricted inho- 
mogeneous perturbations of the form given in Eq. (3.19)). This would mean that 
a macroscopic perturbed FLRW metric with large-scale inhomogeneous terms is 
not compatible with the assumptions that have gone into the theory of MG. In 
particular, MG may only be applicable on the largest scales. 

3.4.1 Discussion 

The first of these possibilities would appear to suggest that perturbed FLRW solutions 
of MG behave in the same way as perturbed FLRW solutions of Einstein's equations 
(to the lowest order of approximation). The MG equations deal with the effects of 
averaging on macroscopic scales, and this is consequently a statement about the be- 
haviour of space-time at, or above, the averaging scale (on smaller scales one would 
expect to have to deal with the linearly perturbed Einstein's equations directly, as 
usual). This result assumes it is appropriate to apply MG in the way that we have 
done in this paper, and that new FLRW solutions to the macroscopic field equations 
that exhibit different behaviour do not remain to be found. If this is the case, then it is 
suggested that on macroscopic scales the FLRW solutions to Einstein's equations are 
a valid description of the averaged geometry (possibly with an 0(e 1 ^ 2 ) renormalization 
of the spatial curvature). This could be regarded as a structural stability result for the 
FLRW solutions of Einstein's equations. 

The second possibility would appear to restrict the domain of applicability of 
MG to macroscopic geometries in which the inhomogeneous perturbations are severely 
restricted, or to situations in which they are spatially homogeneous and isotropic (in 
the simple way that the theory has been applied here, at least). This still leaves 
the possibility of phenomena such as a time-dependent but spatially homogeneous 
anisotropic shear or tilt when describing the universe on macroscopic scales. Such 
effects would appear to be consistent with the anisotropics that arise from averaging 
the flow of matter fields in studies of cosmological backreaction using perturbation 
theory (to lowest order in the shear anisotropy that characterizes the anisotropy of the 
flow) [14], as well as those that arise from a detailed analysis of the observed large 
scale anisotropy of the Hubble flow relative to the CMB rest frame [15]. 

It is clearly on large scales that the effects of averaging are most important, and 
MG can potentially deal with this if perturbations to the macroscopic FLRW geometry 
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are spatially homogeneous. Again, on scales below the size of the averaging domains 
one would expect to be able to use the linearly perturbed Einstein's equations directly. 
What we have found here, however, perhaps suggests that linear perturbation theory 
is not a valid approach on macroscopic scales when A = 0(1). 



4 Distance Measures in the Macroscopic Universe 

Measures of distance are crucial to observational cosmology, as they are often directly 
linked to astrophysical observables such as supernovae and the CMB. Here we will make 
a first step towards calculating angular diameter distances and luminosity distances in 
the FLRW solutions of MG. 

4.1 The Macroscopic Optical Equations 

As is usual in cosmology, we will make use of the geometric optics approximation. The 
behaviour of a bundle of null geodesies in a general space-time is then given by the 
Sachs optical equations, which read [16] 

^+~6 2 -u 2 + a*a= l -r ap k^ (4.1) 

^ + 2uo~6 = (4.2) 

^ + 2a0 = -c a ^ 8 (tT^(t*yk s , (4.3) 

where a is the complex shear scalar, u is the vorticity scalar, and 9 is the expansion 
scalar. The 4-vector k a is tangent to the null curves, and the t a are complex screen 
vectors that are null, orthogonal to k a , with unit magnitude, and that are parallel 
transported along the curves. The variable A denotes an affine parameter used to 
measure position along the null curves. 

In the usual application of Eqs. (4.1)-(4.3) it is common to take r a p = R a $ and 
c a /3 7< 5 = C a/ 3 7< 5, where R a p and C a p^s ar e the Ricci and Weyl tensors of the microscopic 
space-time. Here we are interested in solving Eqs. (4.1)- (4.3) using the average values 
of these quantities. We will therefore take r a/3 = (R a p) and c Q( g 7 5 = (C a ^s)- By 
contracting with the 4- vectors k a and t a , that have been calculated in the macroscopic 
geometry, we can then determine estimates for the driving terms in the evolution 
equations for the expansion and shear of the null congruence (4.1) and (4.3). These, 
in turn, can be used to infer the behaviour of typical values of 9 and a when making 
observations over macroscopic scales in cosmology. 

Using results from [1], we find that 

(R a p) - M aP = Q a p, (4.4) 
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and that 



9 ar (C a ^s) - W^ y5 (4.5) 

= Q T p~t5 - ^ T [ 7 ( 5 fT «5]/3 CT + 9 aT 9li[yQs]a + g 9° P Q a p^ [^9 &}P + ^ f$[^ \S\ 

^ <5[/37 a] ^ 6\p 7tr] ^ a[/3<5 a] ^ 7 [/3 ^ 5[a7/3<r]y ~r ^ 8[a/3ja]9 

1 I 

ycr ycr -olt ryt up T yt pa T 

~ L ~fla8/3a]9 - * i[a06&\9 ~ g Z a [p e]9 ydS/3 ~ 7^ a [ p e ]9 ^95/3 
11111 

± ry^P T - x yep T - ± rvep —f . ± rytp - T i 7« 0"P _ r _ 

+ g Z [pj € ]96p + [ p le] 9S/3 + [p ^ e ]5 7 + 3^ [f>Sle]9 7 + 3^ a[p e]# 5^7/3 

1 1111 

, ye per — f — - 1 - yep r _ - 1 - yep r - x yep -r yep —r 

+ 3^ cr[p e]^ ifl'T/S - [ p 5 e] #7/3 - [ p 5 e] fi , 7/3 - [p^\9 6~ 3 [pj§e\9 5' 

where W T g s is the macroscopic Weyl tensor, constructed from M T /3lS and g a p in the 
usual way. If the macroscopic geometry is FLRW, then we automatically have W a ^s — 
0, so that the average of the Weyl tensor of the microscopic geometry is given only by 
the RHS of Eq. (4.5). The raising and lowering of indices in this equation are done 
using g a/3 and g a p. 

4.2 Optics in Type I Solutions 

Without loss of generality, we choose our coordinate system such that the light rays 
we consider are propagating in the x-direction. We also take the tangent vectors in 
the macroscopic geometry, k^, to be past directed. The null and geodesic conditions, 
k^k^ = and k^k". = 0, then give 

fcM=l(_l, 1,0,0), (4.6) 
cr 

where we have taken the macroscopic geometry to be a spatially flat FLRW solution of 
the macroscopic field equations, and we have used the coordinate system {r],x,y,z}, 
where r) is conformal time (as used in Section 3). A suitable choice of complex screen 
vectors is then 

f* = -L(0,0,l,-i), (4.7) 
V2a 

which are unique up to a rotation in the y-z plane. 

Let us now use these expressions, together with Eqs. (4.4) and (4.5), to calculate 
the driving term on the right-hand side of Eqs. (4.1) and (4.3). These are 



a' 2 a" 



a 2 a 



(4.8) 



and 

-(C al3 ,s)(tT^(tyk 5 = 0, (4.9) 

where, as before, a prime denotes differentiation with respect to rj. That is, the shear 
scalar has no driving term, just as in the FRLW solutions of Einstein's equations, and 
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the total effect of the correlation terms is the presence of A/a on the RHS of Eq. 
(4.1). 

Just as in the macroscopic field equations for spatially homogeneous and isotropic 
geometries, we have that A is the only constant that enters into the final equations. 
What is more, we note that A once again takes exactly the same role as a spatial 
curvature term. That is, if we were to calculate the right-hand side of Eqs. (4.1) 
and (4.3) in a space-time whose microscopic geometry was given by spatially curved 
FLRW, then we would arrive at exactly Eqs. (4.8) and (4.9), with spatial curvature 
k = A 2 . Angular diameter distances and luminosity distances calculated in this way are 
therefore identical to those obtained in a spatially curved FLRW solution of Einstein's 
equations, even though the macroscopic geometry (as found in [3] and [6]) is that of a 
spatially flat FLRW universe. 

4.3 Optics in Type II Solutions 

Let us now consider the RHS of Eqs. (4.1) and (4.3) in the more general Type II 
solutions given in Section 2.2. This time we cannot choose light rays that propagate 
in the x-direction without losing generality. We therefore consider rays with tangent 
vector 

= — (—1, cos 9 cos^, cos 9 smip, — sin 9) , (4-10) 

cr 

which is simply a spatial rotation of the 4- vector k 11 given in Eq. (4.6). Here 9 denotes 
rotation about the y-axis, and ip denotes rotation about the z-axis. Applying this same 
rotation to the screen vectors from Eq. (4.7) gives 

t M = (0, — sinz/) — i cos tp sin 9, cosip — i sinip sin 9, —i cos tp cos 9) . (4.11) 
V2a 



Using these expressions we then find that 



-(R a B)k a k p = 4 

2 \ api fl4 



2— V A cos 9 cos(2^J 

a 2 a 2 



(4.12) 



and 



- (C a ^ s }(t*) a ^(t*yk s = cos 2 g [cos(2^) + zsin(2^) sing] . (4.13) 

These two equations are identical to those of the Type I solutions, given in Eqs. (4.8) 
and (4.9), except for the presence of the terms involving ^{x^y). 

4.3.1 Discussion 

These calculations allow us to consider the effect of inhomogeneities on geometric 
optics, and hence cosmological observations, within the framework of MG. The effects, 
which can in principle depend on both x and y, depend on the direction of the null 
vector k^ and can contribute possible effects both along the null ray and transverse 
to it. Note that when the null ray is entirely orthogonal to the xy— plane, so that 
9 = ±7r/2, there is no additional affect. 
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If either 9 = or ip = then the extra terms involving T on the right-hand side 
of Eqs. (4.12) and (4.13) are identical. This means that if expansion is generated from 
a non-zero J 7 , then the term containing J 7 in Eq. (4.13) will also necessarily drive the 
shear scalar. If ifj = 7r/4, however, then the extra term in Eq. (4.12) vanishes, while 
the imaginary part of the extra term in Eq. (4.13) is non-zero. This means that it is 
also possible to have a driving term in the evolution equation from the shear scalar, 
while not having one in the evolution equation for the expansion scalar. In either case, 
the presence of the terms involving T<i{x y y) in these equations breaks the degeneracy 
between the contribution from the correlation tensor and the inclusion of a spatial 
curvature term. 

As we can choose ^{x, y) to be any function of x, y we like, the amount of freedom 
in the optical properties of these solutions is significant. For example, if 9 = ip = 0, 
then for J z 2 {x,y) > the effect of the extra term is to make distant objects appear 
brighter. This will be due to both the direct consequences of J-~2{x,y) appearing on 
the RHS of Eq. (4.1), as well as indirectly, through its contribution to the evolution 
of the shear in Eq. (4.3). The former of these acts in the same way as Ricci focussing, 
while the latter focusses light rays through the appearance of a*a in Eq. (4.1). 

If 9 = ip = and J-2(x,y) < 0, however, the situation will be somewhat more 
complicated, with the extra term on the RHS of Eq. (4.12) causing de- focusing of 
the light rays, while the contribution from the extra term on the right-hand side of 
Eq. (4.13) will still cause a* a to grow, and hence will cause focusing. It can then be 
the case the de-focusing occurs at smaller distances (before the shear has had time to 
accumulate), while focusing will occur at larger distances (after the shear scalar has had 
sufficient time to grow large, as is generically expected to happen at large distances in 
inhomogeneous space-times). For Ti{x y y) < we can therefore have relatively nearby 
objects appearing dimmer, while more distant objects appear brighter. 

The exact optical properties will of course depend on the precise form of the 
function ^(x, y), as well as the values of 9 and ip. Clearly, if ^{x.y) is allowed to 
have different signs in different regions of the universe, or 9 or ip are allowed to be 
non-zero, then the results could be still more complicated. 

5 Concluding Remarks 

In this paper we have discussed the dynamical evolution equations of the Universe 
on large scales using the theory of Macroscopic Gravity (MG). This theory attempts 
to model the effects of averaging the geometry of space-time and is consequently im- 
portant for the interpretation of cosmological observations. In particular, we have 
investigated spatially homogeneous and isotropic solutions to the MG field equations, 
and presented new exact FLRW cosmological solutions. In these FLRW cosmological 
MG solutions, the macroscopic geometry typically evolves like the FLRW solutions of 
Einstein's equations but with an extra spatial curvature term in the macroscopic equiv- 
alent of the Friedmann equation, even though the macroscopic geometry is spatially 
flat. 
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We then investigated perturbations in the macroscopic geometry around these 
FLRW solutions. We assumed that the scale of these perturbations is larger than 
that of the averaging domains, that the correlation tensor Z a ^ can be perturbed 
as in Eq. (3.6), and that the MG approach is generally valid. We found that the 
macroscopic metric in Eq. (3.1) can only contain spatially homogeneous perturbations 
(or, at most, the very highly restricted inhomogeneous perturbations of the form given 
in Eq. (3.19)). However, as noted earlier, this still leaves the interesting possiblity 
of time-dependent but spatially homogeneous effects such as anisotropic shear or tilt. 
These are permitted on macroscopic scales within MG, and might potentially be of 
interest due to recent studies of cosmological backreaction such as [14] and [15]. 

We then took a first step towards calculating distance measures using the FLRW 
solutions of MG we previously obtained. This was done using the geometric optics 
approximation, and by substituting the averages of the Ricci and Weyl curvature ten- 
sors into the driving terms of the Sachs optical scalar equations. We found that in the 
simplest solutions the correlation tensor acts as a spatial curvature term in the Sachs 
equations, in the same way that it does in the MG field equations. In more general- 
ity, we found that the optical properties depend on terms involving the free function 
^{Xyi/). These terms appears on the right-hand side of Eqs. (4.12) and (4.13)), and 
break the degeneracy between the contribution from the correlation tensor and the 
inclusion of a spatial curvature term. They also depend on the direction of the null 
rays to which is tangent. 

Recent studies on the effect of inhomogeneities on the optical properties of uni- 
verses that are statistically spatially homogeneous and isotropic on large scales have 
suggested that their consequences can be large enough to be of importance for inter- 
preting observational data (see, for example, [17]). The motivation for the current work 
is to start detailed and explicit computations of the effects of inhomogeneities on cos- 
mological observables, such as the growth of large-scale structure and the propagation 
of light rays, within the framework of MG. 

Appendix A: Type II Macroscopic FLRW Solutions 

Here we shall generalize and expand upon the solution to the MG equations found by 
van den Hoogen in [3]. To recap, the assumptions made in [3] are 

1: We macroscopic metric g a/3 is the assumed to obey the splitting rule {gapl 1 5e) = 

9afTpSe- 

2: The line-element of the macroscopic geometry can be written 

ds 2 = g a /sdx a dx 13 = a 2 {rj)[—drf + dx 2 + dy 2 + dz 2 \ 

and there exists a time-like unit vector field orthogonal to the spatial hyper- 
surfaces of homogeneity and isotropy, u a = ^-(—1, 0, 0, 0). 

3: The averaged microscopic energy-momentum tensor can be modeled macroscop- 
ically as a perfect fluid with energy density p and pressure p. 
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4: The electric part of the connection correlation tensor is zero, i.e. Z a ^ iycr u a = 0, 
and, finally, 

5: The connection correlation tensor Z a ^ va and the affine deformation tensor 
A a pj are invariant under the same of Killing vectors as the macroscopic 
metric. 

We note that if one interprets the contribution of the connection correlations in the 
MG equations to be an effective fluid with isotropic pressure p gra v and energy density 
Pgrav, then assumption 4 implies p grav + 3p grav = 0. In what follows we shall assume 
only that the first 4 assumptions listed above hold true, while relaxing assumption 5. 

In general, without taking into account the cyclic constraint, the connection cor- 
relation tensor has 720 independent components. Assumption 4 immediately sets 360 
of them to zero. We shall label the remaining 360 variables as fi(x,y, z,rj). This 
assumption is critical to solving the quadratic constraint (1.14). The equi- affine con- 
straint and the cyclic identity equations (1.11) and (1.10) then yield 239 constraints, 
leaving 121 independent variables remaining. The integrability condition (1.13) then 
yields an additional 52 constraints, leaving 69 independent variables remaining. It is 
interesting to note at this point that the polarization tensor Q a p fll/ contains only 18 
of these 69 independent variables, and Q a /3fJja = its trace, contains only 6. The 
differential equations (1.12) now constrain the remaining independent 69 variables of 
Z a p ^ to be functions of x, y and z only. The remaining equations yield a set of 40 
linearly independent differential equations in the spatial variables. 

We shall now replace assumption 5 with the following set of less restrictive as- 
sumptions: 

5a: Only SitGT 9 ™^ = -{Z\ V1 - ^Q^g^, is invariant under the £ 6 group of 
motions. 

5b: The 40 independent variables A a p are all assumed to be independent of z. Eq. 
(1.15) then implies that the 18 independent variables that determine Q a ^ (a 
subset of the fi(x, y, z)) are also independent of z. 
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5c: We assume that: 



A 1 

UU 


= o, 


i G {1,2,3} 




= o, 


i G {1,2,3} 




= o, 


G {[1,2], [2, 3], [3,1]} 


A { 


= 0, 


[ij] G {[1,2], [2,1], [2, 3], [3, 2], [3,1], [3,1]} 


71 33 


= 0, 


ze {1,2} 




= 0, 


ie{l,2} 




= 0, 


ie{l,2} 


A 1 


= 0, 


ze{i,2} 


A i 

^jk 


= 0, 


[i, j, fc] G {[1,2, 3], [2, 3,1], [3, 1,2]} 


A 1 


= AJ 0ji 


{[1,2], [2, 3], [3,1]} 




= A \i, 


M e {[1,2], [2,1]} 



Assumption 5a yields 5 algebraic constraints and 3 differential constraints, so that 
there are now 64 independent variables in Z° ' p^ V(J - We note that Q a ' p^ u contains only 
13 of these 64 independent variables, and Qp^ contains only 6. The number of linear 
independent differential equations reduces to 37. One may note that a solution to 
these 37 differential equations is to let each of the remaining independent variables in 
Z a pJ l V(T be constant. 

Assumption 5c yields 30 constraints on the A a g , leaving 10 independent variables 
in A a p . Assuming 5b and 5c together, and solving Eqns. (1.15) and (1.16) simulta- 
neously, reveals an additional 11 constraints. This leaves 53 independent variables in 
Z a o J 1 ^ of which only 2 are found in Q a p^ u and Qp^. We label these two "special" 
independent variables as Ti(x,y) and T2(x,y). These assumptions also decrease the 
number of differential equations determining Z a ^ ya by 5. Surprisingly, the two func- 
tions J-"i(x, y) and ^{x, y) are not restricted or determined at all by what is now a set 
of 32 differential equations for the other 51 variables fi(x,y,z) in Z a ^ V(y . One can 
now solve both Eq. (1.15) and Eq. (1.16) to find four families of solutions parameter- 
ized by the ordered pair (a 1} a 2 ) = (1/2, 0), (1/2, 1/2), (1, 0), (1, 1/2). The values for 
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A a p 1 in these families are 



AX 

XX 


= (1- 


a 2 )Bi 


y) 


A y 

XX 


= (1 - 


a 1 )B 2 \ 


'x, v) 


AX 

xy 


= (1" 


ai)B 2 \ 


fay) 


A^ m , 
xy 


= a 2 Bi 


(x, y) 




AX 

A yy 


= a 2 Bi 


fay) 




AV yy 


= OL\B 2 


fay) 




A\ z 


= B 3 






A\- 


= A 






A i 

A i V ■ 


= -A 






4° 

A 00 


= -A, 







which depend on two constants, A and B 3 , and two functions, Bi(x,y) and B 2 (x,y). 
The functions Bi(x,y) and B 2 (x,y) must satisfy the following system of differential 
equations, which are the remaining equations from (1.15), 

(a 2 - ai)^j- = 7i (5.1) 

(a 2 - = ft ( 5 - 2 ) 

dB 1 ,dB 2 , . 

h [Oil ~ 1)-^— = -T2, (5.3) 

ox ay 

where the functions T\{x,y) and T^x^y) are the two "special" functions found in the 
polarization tensor Q a p fiu - If a± ^ ct 2 then these three differential equations can be 
considered as definitions for Fi(x,y) and T 2 (x,y), given any two functions B\(x,y) 
and B 2 (x,y), that satisfy = If ai = t*2, then F^x^y) = 0, and the third 

equation can be interpreted simply as a definition for ~F 2 (x,y). For completeness, the 
values of the polarization tensor Q a l3 ^ v are 

Q xyx Q V yxy J~l(x, y) 

Q\ yx = -A 2 + F 2 fay) 

Q x yxy = A 2 + T 2 (x,y) 

O x = O v = O z = O z = A 2 

^ zzx V zzy xxz V • 

We will not explicitly present the correlation tensor Z a n J M va here, as it has 1584 non- 
trivial components in the 54 variables fi(x,y,z), F\{x,y), F 2 {x,y) and A. 
It should be noted that these solutions give 



SvrGT^ = — — diag(3, 1,1,1) 



A 2 

a{r]f 

which is exactly the same as the Type I solutions. The contibution of the correlation 
tensor to the MG field equations is therefore the same for both the Type I and Type II 



- 20 - 



solutions, despite the latter having considerably less symmetry in Z a ^ V(J and A a / g 7 . 
Finally, one should also note that it may be possible to find even more solutions if 
assumptions 5b and 5c are relaxed. 
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